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ABSTRACT 



The perfomance of a phased array antenna depends upon the phase and 
amplitude distribution acax)ss the ^)erture. A continuous linear phase 
distribution is required for an efficient focussed beam, but practical 
phase shifters are digital devices and can only provide an approximation 
to a linear phase distribution. Consequently seme type of roundoff 
criteria most be established. The method of roundoff affects the radiation 
pattern characteristics such as beam location and sidelobe level. Accurate 
target tracking requires that a radar have a small beam pointing error. 
low sidelobes are also desirable to prevent jaimdng and the illumination 
of clutter. Therefore the goal was to select a roundoff criterion that 
provides a phase distribution across the ^>erture to minimize both the 
beam shift and the sidelobe levels, vhile simultaneously maximizing the 
gain. The methods examined are referred to regular roundoff, wei^ted 
randem roundoff, running sum roundoff and symmetric running sum roundoff. 
The first two are in oemron use, but the third and fourth are new methods 
examined in this paper. It was daronstrated that the latter two have the 
minimum beam shift of the four roundoff methods, without significantly 
degrading the other pattern properties. 
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I. 



INTRODUCTION 



Early radar and communication systems used antenna arrays 
consisting of a combination of separate radiating elements. 
These systems date back to the turn of the century [Ref. 1: p. 
11-2]. The antenna radiation characteristics are determined by 
the position of the radiating elements and the amplitude and 
phase excitation relationships. As the operating frequency of 
radar systems increased (decreased wavelength), array antenna 
radars where replaced by simpler parabolic antennas. With the 
advent of printed circuit technology, electronically 
controlled phase shifters have become more compact and low 
cost. Consequently attention has been directed to array 
antennas again. The aperture excitation may now be modulated 
by controlling the individual elements to give beams that are 
scanned electronically [Ref. 1: p. 11-2]. An electronically 
steered array antenna can track a large number of targets, 
illuminate targets with RF energy to guide missiles toward 
targets, perform a complete hemispherical search with 
automatic target selection, and perform multiple functions in 
communication systems [Ref. 1: p. 11-2]. 

The performance of a phased array antenna depends upon the 
phase and amplitude distribution across the aperture. A 
continuous linear phase distribution is required for an 
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efficient focussed beam, but practical phase shifters are 
digital devices and can only provide an approximation to a 
linear phase distribution. Consequently some type of roundoff 
criteria must be established. However the method of roundoff 
affects the radiation pattern characteristics such as beam 
location and sidelobe level. Accurate target tracking requires 
that a radar have a small beam pointing error. Low sidelobes 
are also desirable to prevent jamming and the illumination of 
clutter. Therefore the goal is to select a roundoff criterion 
that provides a phase distribution across the aperture to 
minimize both the beam shift and the sidelobe levels, while 
simultaneously maximizing the gain. 

This paper examines the performance degradation of a linear 
array due to the quantization of phase, which is present in 
all practical electronically scanned arrays. Analysis was 
restricted to linear arrays so that the simulations were 
computationally fast. The same techniques can be extended to 
two and three dimensional arrays. The performance of a 
monopulse antenna were simulated on the computer for a linear 
array using the Bayliss distribution. A Bayliss distribution 
is a low sidelobe difference beam distribution commonly used 
in radar. An ideal difference beam has a perfect null at the 
scan angle for a continuous linear phase distribution. When 
phase and amplitude errors are present, a perfect null is not 
possible. In practice a good null is about 35 dB below the sum 
beam peak [Ref. 1: p. 21-22]. 
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The size of the arrays studied was be limited to 26, 50, 
and 7 6 radiating elements using 3, 4, and 5 bit phase 
shifters. A phase bit is a measure of how accurately a phase 
shifter can represent a linear phase by a step approximation; 
the more bits, the closer the phase steps. Although a large 
number of bits reproduce a linear phase better, they add cost 
and complexity to the antenna. Thus there is a tradeoff 
between the bitsize and antenna performance. In addition to 
the phase bitsize, the method of roundoff also affects the 
antenna performance. Four different methods of phase roundoff 
were examined to access their impact on antenna performance. 
The methods examined are referred to as regular roundoff, 
weighted random roundoff, running sum roundoff and symmetric 
running sum roundoff. The first two are in common use, but the 
third and fourth are new methods examined in this paper. It 
was demonstrated that the latter two will have the minimum 
beam shift of the four roundoff methods, without significantly 
degrading the other pattern properties. The antenna system 
parameters evaluated include the beam pointing error, maximum 
sidelobe level, and the null depth of the difference beam. 

This paper is divided into the following sections: basic 
antenna array theory, phase shifters, phase quantization 
algorithms, and analysis of the data. 
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II. PHASE SHIFTERS: AN OVERVIEW 



The evolution of sophisticated phased array antennas is due 
to two factors: the development of electronically controlled 
phase shifters and the integration of computers into the 
system. Electronically controlled phase shifters allow the 
beam to be redirected in less than a few microseconds and the 
computer can process a large amount of data. The principle 
advantage of these two factors is the minimum amount of time 
spent moving a beam between targets and making decisions about 
detection and tracking [Ref. 1: p. 11-2]. The performance of 
a phased array antenna is highly dependent upon its phase 
shifters. An N bit phase shifter will have a most significant 
bit of n radians (180 degrees) and a least significant bit of 
2n/2^ radians. For best performance, a large number of bits 
should be used in the phase shifters, while for minimum cost 
and complexity, a small number of bits is desired. Since most 
feed methods require a phase shifter at every element of the 
array, the performance, size and cost of the antenna system 
are strongly affected by the phase shifter [Ref. 1: 12-3,4]. 

Several parameters of phase shifter design impact the 
performance of the phased array radar. The issues of concern 
are insertion loss, switching time, drive power, phase error, 
and physical size. The insertion loss should be as low as 
possible to reduce the loss of power during transmission and 
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lower the signal-to-noise ratio at the receiver. A low 
insertion loss also reduces phase shifter heating due to ohmic 
losses. If amplifiers are used between the phase shifter and 
the radiating element, the power capacity and insertion loss 
become insignificant. The phase shifter switching time should 
be as short as possible. A long switching time increases the 
minimum radar range when non-reciprocal phase shifters are 
used and when a burst of pulses are transmitted in different 
directions [Ref. 1: p. 12-3,4]. 

The drive power should be as small as possible. A large 
amount of drive power generates heat and may require power 
supplies that are too large for the given system. They also 
require expensive driver circuit components [Ref. 1: p. 12-3]. 
The phase error should also be as small as possible. It should 
not reduce the antenna gain substantially on transmitting or 
raise sidelobes on receiving. One cause of phase error is the 
size of the least significant bit of digital phase shifters. 
Others are due to the tolerances in the phase shifters and 
driver. The maximum acceptable physical size is limited by the 
element separation (usually A/2 by A/2) in order to be 
packaged behind each radiating element [Ref. 1: p. 12-3,4]. 

Today there are two types of phase shifters used in radar 
systems: ferrite and diode. Ferrite phase shifters are analog 
devices, although in practice they are digitally controlled. 
The first ferrite phase shifter employed in an operational 
phased array was the Reggia-Spencer phase shifter. The 
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development of a latching toroidal ferrite phase shifter, 
which employs a transverse magnetizing field, made significant 
improvements . The driving power was reduced and the driver 
circuitry was simplified. Also, the switching time was reduced 
from about a millisecond to a few microseconds. Ferrite phase 
shifters are impractical in radar systems operating at L-band 
(1 to 2 Ghz ) frequencies and below since the magnetizing field 
is operating below resonance. This is due to available 
materials, which have a lower limit of saturation moment of 
approximately 200 gauss. At L band, a large saturation moment 
produces peak power limiting at a low power level, in a 
reasonably compact phase shifter. In addition, a low moment 
material is very temperature sensitive [Ref. 1: p. 12-4,5]. 

Diode phase shifters were developed with the recognition 
of the importance of the pn diode in radar systems . As 
fabrication techniques improved the diode breakdown voltage, 
cutoff frequencies, and reliability increased. Diode phase 
shifters are temperature insensitive over a wide range and do 
not require air or water cooling, as do the ferrite phase 
shifters. Use of diode phase shifters in radar systems is 
practical at all operating frequencies. At lower frequencies 
(below L band) diode phase shifters perform better since the 
semiconductor loss decreases as frequency decreases. This 
decrease can be traded for high peak power performance by 
adjusting the impedance of the transmission line on which the 
diodes are mounted [Ref. 1: p. 12-4,5]. 
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At frequencies at or above L band, the selection of the 
phase shifter type is more difficult, because of high diode 
phase shifter losses and the difficulty in maintaining 
manufacturing tolerances. Table 1 summarizes the performance 
of the phase shifters at these frequencies [Ref. 1; p. 12-5]. 



TABLE 1: SUMMARY OF PERFORMANCE OF TYPICAL FERRITE AND DIODE 
PHASE SHIFTERS 





FERRI! 

PHASE 


PE (TOROIDAL) ::r- 
SHIFTER : 


DIODE PHASE 
SHIFTER; 


INSERTION LOSS 


1.0 dB 


1. 0-2.0 dB, 
increases with 
the number of 
bits 


SWITCHING TIME 


2. 0-5.0 psec 


50 nsec - 2 psec 


DRIVE POWER 


2 W AT 1,000 pps 


1.0-2. 5 W 


TRANSMITTED POWER 


75 kW PEAK, 
400 W AVERAGE 


10 kW PEAK, 
200 W AVERAGE 
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III. ANTENNA ARRAY THEORY 



In the following discussion an antenna array refers to a 
group of similar elements arranged in an arbitrary geometric 
configuration, and excited with a prescribed phase and 
amplitude relation to give a desired radiation pattern. The 
radiation characteristics most commonly examined are the 
direction and width of the main beam, sidelobe levels, and 
null depth or main lobe level. In this chapter the 
fundamental theory and characteristics of linear arrays will 
be examined. 

The total radiation pattern is the sum of the contributions 
from the individual elements. For a linear array having N 
radiators, with element n excited by amplitude, a^^, and phase, 
i|jr^, and inter-element spacing s, the radiation pattern is 



w-i 






ji^nsieiniQ) ■ 






f{0) 



(lii-l) 



n=0 



where f(0) is the radiation pattern of a single element [Ref. 
1: p. 11-9]. A typical linear array is shown in Figure 1. The 
term includes the inter-element phase to scan the beam as 
well as any phase errors that may be present. For convenience, 
the scanning phase of 2ns s in (0^^) /A. radians per element will 
be separated out so that is entirely due to the phase error 
[Ref 1. : p. 11-9], 
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E{d) =f{d)e'^~^ 



2)t [w-D 

’T“1— 



s(sln(6) -einOj) 




ji[-^/xs(sin{0) -sin ( 00 ) 



(III-2) 



The leading exponential references the phase to the center of 
the array and will be dropped since it represents a phase 
shift conunon to all elements and does not affect the antenna 
pattern . 



Figure 1 . Linear Array with N Radiators 
Uniformly Spaced by a Distance s 

If there are no phase errors and all amplitudes are equal 
(i|r^=0 and a^= 1), then equation 1 1 1-2 can be reduced to a 
closed form expression, [Ref 1; p. 11-9] 




sin (sin (0) +sin (Og) ) ] 



£:(e)=f(0) 

i\^sin[7i-^ (sin(0) +sin(0g) ) ] 



(III-3) 
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Equations I I 1-2 and I I 1-3 are the fundamental response 
equations for a scanned array system. As a typical example, 
Figure 2 shows the radiation pattern calculated from equation 




III-2 , with 50 elements and s = ,3X . In this case all are 
zero and the a^ are chosen to give a 25 dB Taylor 
distribution. The radiation pattern will have only one major 
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lobe and grating lobe maxima will not occur for values of -rr/2 
< 0 < n/2 as long as [Ref. 1: p. 11-11] 

(sin(0) -sin(0J )<±7i . (III-4) 

The beam of an antenna points in the direction normal to the 
phase front. In phased arrays, the phase front is adjusted by 
controlling the phase of each radiating element, as shown in 
Figure 3 . 




A. RADIATING ELEMENTS 

The most common type of radiating elements used in array 
antennas are slots, dipoles, horns, and more recently 
conductive patches. Since the size of the radiating element is 
constrained to fit in the array geometry, the area of the 
elements are limited to approximately The predominate 

electrical factors in choosing the radiating element in an 
array environment is the impedance and the pattern of a 
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radiator. Another factor that influences the type of element 
is the feeding method. The radiating element should be chosen 
so that the antenna aperture and feed system can be packaged 
to satisfy the size and weight reguirements of the antenna. 
For example, if the feed system contains stripline 
components, a stripline dipole would be the optimal choice. On 
the other hand, if waveguide feed devices are used, open-ended 
waveguides or small horn antennas are more convenient [Ref 1: 
p. 11-5]. 

Most well matched antenna elements have element patterns 
that can be approximated by f(0) = cos(0), at least for angles 
near broadside (normal to the antenna aperture) . In this 
analysis only the array factor is considered; the element 
pattern will be neglected (f(0) = 1). This simplification will 
not significantly affect the results because the scan angles 
of interest are limited to 0^^ = ± 60°, and comparisons are 
referenced to an error-free array also with f(0) = 1. 
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IV. PHASE QUANTIZATION ALGORITHMS 

In order to steer a focused beam to a certain direction, 
a linear phase must be set up across the aperture according to 
equation I I 1-2. However, most phase shifters are digital 
devices and the allowable phase values are limited by the bit 
size. Therefore, in most instances the phase at each element 
must be quantized. For a n bit phase shifter the bit size is 
determined by 

Bit 5ize = -^^ . (IV-1) 

2 ” 

Each allowable phase shifter value is referred to as a phase 
state. The phase states are given by 

{i-1) ‘Bit Size i=l , 2 , 3 , . . . . , 2” . (IV-2) 

For example, a 4 bit device has a bit size of 360°/2^=22 . 5° and 
16 phase states. The difference between the desired phase and 
the quantized phase is denoted the quantization error. 

In this paper four different methods of phase quantization 
were used to examine the phased array antenna performance. 
They are regular roundoff, weighted random roundoff, running 
sum roundoff and symmetric running sum roundoff. Each of these 
methods of phase quantization is discussed in detail in the 
following sections. 
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To examine the impact of quantization error on the 



difference beam radiation pattern consider equation III-2: 



^ j (3pS n (8in(6) -sin(ea) ) 

/ 3=0 



(IV-3) 



where a^ is the Bayliss coefficient and is the quantization 
error at the nth element. Equation IV-3 assumes a line source 
and that the only errors present are those due to the phase 
quantization. Therefore at any scan angle 



N-l 

E{d^)=Y^a„e^'*° . (iv-4) 

n=0 

For a difference beam the coefficients a^, will be positive for 
one half of the elements (1 < n < N/2 , N even) and negative 
for the other half (N/2 + 1 < n < N) . Thus, if a null is to 
occur at 0^^ then |E(0 q) 1 -+ 0, [Ref. 2: p. 2] and the weighted 
sum of the quantization errors must be zero. 

A . REGULAR ROUNDOFF 

Regular roundoff is the most straightforward of the four 
methods considered. It involves calculating the desired phase 
at each radiating element and then rounding the phase up or 
down to the closest allowable state of the phase shifter. This 
method of phase estimation results in a periodic quantization 
error, the period of which varies with the bitsize and scan 
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angle. A typical desired versus quantized phase angle for 
regular roundoff is shown in Figure 4. This type of phase 
error causes large sidelobes and potentially large beam 
pointing errors. 
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Figure 4. Desired (o) and Regular Roundoff (*) Phase Angle 
Distribution for a 26 Element Array, 3 Bit Phase Shifter 

B. WEIGHTED RANDOM ROUNDOFF 



This method rounds up or down randomly such that the 
expected value of the phase at any element is the desired 
value. This destroys the periodic quantization error, but also 
raises the average sidelobe level. 
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This method is implemented by choosing a random number 
uniformly distributed between 0 and 1 and assigning it to an 
element of the array. The algorithm combines the random number 
and the desired phase to determine the phase state at the 
given radiating element as follows. At each element the 
algorithm subtracts the closest phase state from the desired 
phase (which gives the regular roundoff quantization error) 
and divides the difference by the bitsize. The algorithm then 
adds this to the random number chosen for this particular 
element. If the sum is greater than one, the phase at the 
radiating element is rounded up, otherwise the phase is 
rounded down. Using this method, the expected value of the 
quantization error tend towards zero as the number of the 
array elements increases. A typical desired versus weighted 
random roundoff algorithm phase angle, i|x, is shown in Figure 
5 . 



C. RUNNING SUM ROUNDOFF 

This method chooses the phase state at the nth radiating 
element as the one which keeps the magnitude of the cumulative 
sum of the quantization errors the smallest. This minimizes 
the cumulative quantization error at any point along the 
array. For instance, at element n the algorithm compares the 
cumulative quantization error that results from either 
rounding up or down at the element. The algorithm then chooses 
the phase state that minimizes this total quantization error. 
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This condition will give |E(0^)| 0 as the number of elements 

is increased, A typical desired versus running sum roundoff 
algorithm phase angle, ij;, is shown in Figure 6. 

D, SYMMETRIC RUNNING SUM ROUNDOFF 

This method is related to the running sum roundoff just 
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Figure 5. Desired (o) and Weighted Random Roundoff (*) Phase 
Angle Distribution for a 26 Element Array, 3 Bit Phase Shifter 

described, but takes advantage of the phase symmetry across 

the array. If the phase reference is defined at the center of 



17 



the array, then elements located at equal distances from the 
center will have equal but opposite phase values. For one half 
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Figure 6. Desired (o) and Running Sum Roundoff Phase Angle 
Distribution for a 26 Element Array, 3 Bit Phase Shifter 



of the array, the running sum roundoff algorithm can be used 
to determine the quantized phase. The phase state for the 
corresponding elements on the second half of the array will be 
set so that the quantization error is equal and opposite. 

The rational for this approach is based upon equation IV- 
4, which for a sum beam with a symmetric amplitude 
distribution becomes. 
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N 

2 



=E^n( 






(IV-5) 



For small errors, 



N 

2 

E(0o)*Ean(2+jtn + j'l^*1-T,) • (IV-6) 

n=1 

Thus choosing will force the error in each term of 

the sum to zero. In practice this is not possible because the 
cosine term from the exponential, which leads to the 2, is 
significantly in error. Therefore the peak of the beam is 
somewhat less than S a^. This method is included in the study 
to determine whether there is a symmetric property that can be 
exploited. A typical desired versus symmetric running sum 
roundoff algorithm phase angle, ijf, is shown in Figure 7. 

E. QUANTIZATION ERROR AND NULL DEPTH AND LOCATION 

In the case of a difference beam, if the roundoff 
algorithm could set up a linear phase across the aperture at 
a scan angle, then there would be a perfect null at that angle 
(i.e. zero pointing error and zero gain). Since the bitsize of 
the phase shifter is limited, applying a roundoff algorithm to 
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an array fills the null and introduces a pointing error, as 
shown in Figure 8. The phase error causes the null to be 
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Figure 7. Desired (o) and Syirunetric Running Sum Roundoff Phase 
Angle Distribution for a 26 Element Array, 3 Bit Phase Shifter 



rounded instead of a sharp notch, and moves the minimum value 
away from 0^. Both the location and depth of the null are 
important. If a null is sharp but located at the wrong angle, 
a tracking error will be introduced in a radar system. If the 
null is too shallow (rounded), then the radar will not be able 
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to position the minimum on the target, and again a tracking 
error will be introduced. 
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Figure 8. Perfect Bayliss Difference Beam Null versus Roundoff 
Algorithm Null 
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V. ANALYSIS OF THE DATA 



In this chapter, data from the different phased array 
antenna configurations will be examined. Comparisons will 
examine the effects on the null depth, sidelobe levels, and 
pointing error as a function of the number of elements, 
roundoff algorithm, and the number of bits. Only one of the 
these will be changed in any given comparison. 

A. COMPARISON BY ROUNDOFF ALGORITHMS 

In this section, typical data is presented for an array 
of 50 elements using the four roundoff algorithms with 3 
bit phase shifters. Figures 9 through 12 illustrate the range 
of sidelobes, null depth and null location that results from 
quantization errors. The statistics from a large number of 
calculated patterns at 61 scan angles (0® to 60® in 1® steps) 
are presented in tables 2 through 7 . 

The running sum and symmetric running sum roundoff 
algorithms had the minimum pointing error, in most cases a 
factor of 10 smaller than the other two algorithms. All four 
algorithms were comparable as far as sidelobe levels are 
concerned . 
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Figure 9. 50 Element Array Bayliss Difference Beam Radiation 
Pattern for a 3 Bit Phase Shifter using Regular Roundoff 
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Figure 10. 50 Element Array Bayliss Difference Beam for a 3 
Bit Phase Shifter using Weighted Random Roundoff 
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Figure 11. 50 Element Array Bayliss Difference Beam for a 3 
Bit Phase Shifter using Running Sum Roundoff 
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Figure 12. 50 Element Array Bayliss Difference Beam for a 3 
Bit Phase Shifter using Symmetric Running Sum Roundoff 
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TABLE 2. COMPARISON OF 26 ELEMENT ARRAYS WITH 25 dB SIDELOBE 



LEVEL USING DIFFERENT ALGORITHMS 



Number of 
EleBoeats . 
/Roundoff 
Algorithm 


# Of 
Bits ; 


Pointing Error 
Mean / Variance 
< milli Degree) 


Null 

Dej^th 

Mean 

<dB) 


Side lobe 
Level Mean/ 
Variance (dB) 


26/ REG 


3 


23.10 / 66.84 


-66.39 


-17.24 / 7.87 


26/ WGT 


3 


-10.10 / 123.72 


-44.03 


-16.00 / 3.80 


26/ RUN 


3 


- 4.67 / 3.97 


-63.00 


-17.66 / 8.69 


26/ SYM 


3 


- 1.64 / 3.83 


-88.64 


-15.53 / 6.48 


26/ REG 


4 


7.85 / 13.59 


-77.49 


-21.03 / 3.76 


26/ WGT 


4 


-44.44 / 34.40 


-52.75 


-20.10 / 2.49 


26/ RUN 


4 


0.90 / 0.91 


-75.75 


-20.95 / 3.09 


26/ SYM 


4 


- 2.39 / 1.11 


-88.27 


-19.09 / 4.56 


26/ REG 


5 


- 9.93 / 5.11 


-82.72 


-23.56 / 1.10 


26/ WGT 


5 


-13.80 / 8.72 


-66.90 


-22.80 / 1.03 


26/ RUN 


5 


- 0.31 / 0.27 


-82.22 


-23.20 / 0.80 


26/ SYM 


5 


1.70 / 0.34 


-87.02 


-22.39 / 2.42 



TABLE 3. COMPARISON OF 26 ELEMENT ARRAYS WITH 40 dB SIDELOBE 



LEVEL USING DIFFERENT ALGORITHMS 



Number of ^ 
El^tents 
/Roundoff 
Algorithm 


# of ^ 
Bits ; 


Pointing Error 
Mean / Variance 
(mini Degree) 


Hull 

Dejpth 

Mean 


Side 1 obel^l^iM 


Level Mean / 
Variance (dB) - 


26/ REG 


4 


6.74 / 10.760 


-75.94 


-20.95 / 5.69 


26/ WGT 


4 


- 8.82 / 26.302 


-52.41 


-21.52 / 8.25 


26/ RUN 


4 


- 2.16 / 1.525 


-62.09 


-21.22 / 6.86 


26/ SYM 


4 


- 4.00 / 2.370 


-87.43 


-20.76 / 5.63 


26/ REG 


5 


-10.80 / 4.981 


-83.64 


-28.55 / 7.78 


26/ WGT 


5 


- 8.77 / 11.946 


-66.50 


-26.89 / 3.33 


26/ RUN 


5 


0.64 / 0.376 


-66.56 


-26.54 / 4.15 


26/ SYM 


5 


3.11 / 0.684 


-91.24 


-25.79 / 3.42 
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TABLE 4. COMPARISON OF 50 ELEMENT ARRAYS WITH 25 dB SIDELOBE 
LEVEL ARRAY USING DIFFERENT ALGORITHMS 



Nuinber : 

Sleatenl^s 
/Roundoff i 
Algorii:hm ^ 


# Of ^ 
Bits ; 


Pointing Error 
Mean / Vnrianoe 
(MiiXi Begreo) 


mil 

Depth 

;Kean 

<4B) 


SideXobe 
LeveX Mean / 
Variance 


50/ REG 


3 


- 0.31 / 9.337 


-68.52 


-17.38 / 5.44 


50/ WGT 


3 


-32.85 / 18.455 


-48.25 


-17.56 / 2.38 


50/ RUN 


3 


- 0.34 / 0.191 


-65.69 


-17.75 / 6.01 


50/ SYM 


3 


- 2.69 / 0.191 


-81.46 


-16.18 / 6.48 


50/ REG 


4 


- 5.74 / 5.126 


-76.98 


-21.99 / 2.69 


50/ WGT 


4 


1.82 / 2.982 


-56.64 


-21.76 / 1.50 


50/ RUN 


4 


- 0.85 / 0.052 


-75.30 


-21.88 / 2.53 


50/ SYM 


4 


2.06 / 0.055 


-83.68 


-21.21 / 2.54 


50/ REG 


5 


- 0.24 / 0.589 


-81.49 


-24.24 / 0.78 


50/ WGT 


5 


4.50 / 1.616 


-68.57 


-23.67 / 0.88 


50/ RUN 


5 


1.08 / 0.011 


-80.06 


-24.20 / 0.40 


50/ SYM 


5 


0.60 / 0.017 


-81.18 


-24.14 / 1.05 



TABLE 5. COMPARISON OF 50 ELEMENT ARRAYS WITH 40 dB SIDELOBE 



LEVEL ARRAY USING DIFFERENT ALGORITHMS 



Nuxaber of 
:: EX^tentS 
/Roundoff 
Algorithm : 


# Of 

Bits 


Pointing Error 
Mean / variance 
» <miXli Degree) 


HuXX 

Depth 

Mean 

my 


SideXobe 
LeveX Mean / 
Variance <dB) 


50/ REG 


4 


- 4.02 / 4.210 


-76.49 


-21.01 / 7.77 


50/ WGT 


4 


14.80 / 5.863 


-58.41 


-21.89 / 5.99 


50/ RUN 


4 


- 0.62 / 0.069 


-58.30 


-23.09 / 6.23 


50/ SYM 


4 


3.38 / 0.111 


-84.18 


-22.32 / 6.79 


50/ REG 


5 


0.84 / 0.839 


-82.35 


-29.17 / 5.21 


50/ WGT 


5 


8.54 / 2.078 


-67.45 


-28.83 / 1.91 


50/ RUN 


5 


0.67 / 0.014 


-IQ. 61 


-28.59 / 3.20 


50/ SYM 


5 


0.77 / 0.036 


-83.14 


-27.19 / 5.12 
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TABLE 6. COMPARISON OF 76 ELEMENT ARRAYS WITH 25 dB DESIRED 
SIDELOBE LEVEL USING DIFFERENT ALGORITHMS 



Humber of 
Slements 
/Roundoff 
Algorithm 


# Of : 


Poiutiug Error 
Mean / Variance 
(miiii Degree) 


HuXI 
Depth 
Mean . 
<dB) 


SideXobe 
Level Mean / 
Variance (dB) 


76/ REG 


3 


4.18 / 9.655 


-70.19 


-17.98 / 4.95 


7 6/ WGT 


3 


- 5.28 / 4.117 


-48.50 


-18.66 / 2.50 


76/ RUN 


3 


- 0.30 / 0.042 


-68.41 


-18.15 / 6.99 


76/ SYM 


3 


- 1.15 / 0.050 


-80.79 


-16.98 / 6.24 


76/ REG 


4 


- 2.28 / 0.142 


-74.74 


-22.59 / 3.10 


7 6/ WGT 


4 


- 0.41 / 1.863 


-58.29 


-22.46 / 1.05 


76/ RUN 


4 


- 0.06 / 0.007 


-74.55 


-22.51 / 1.65 


76/ SYM 


4 


- 0.07 / 0.008 


-78.31 


-22.64 / 2.18 


76/ REG 


5 


- 2.69 / 0.112 


-78.05 


-24.42 / 0.60 


76/ WGT 


5 


- 4.98 / 0.466 


-68.49 


-24.04 / 0.86 


76/ RUN 


5 


- 0.05 / 0.002 


-78.78 


-24.22 / 1.91 


76/ SYM 


5 


0.05 / 0.002 


-79.54 


-24.61 / 0.57 



TABLE 7. COMPARISON OF 76 ELEMENT ARRAYS WITH 40 dB SIDELOBE 



LEVEL USING DIFFERENT ALGORITHMS 



^•EX<^nta #■ 
■J^gorithm'®^ 


■I'of " 
'Bits 


: Pointing Error 
Mean / Variance 
(miXii Degree) 


HuXX 

Depth 

Mean 

m) 


SideXobe 
Level Mean./ 
Variance (dB) 


76/ REG 


4 


- 3.66 / 0.213 


-78.14 


-24.13 / 4.80 


7 6/ WGT 


4 


- 1.80 / 1.868 


-58.55 


-25.31 / 3.36 


76/ RUN 


4 


0.43 / 0.007 


-64.40 


-24.19 / 4.77 


76/ SYM 


4 


- 0.18 / 0.017 


-80.44 


-23.66 / 4.39 


76/ REG 


5 


- 2.96 / 0.123 


-79.24 


-29.68 / 2.42 


76/ WGT 


5 


0.36 / 0.571 


-68.22 


-30.18 / 2.10 


76/ RUN 


5 


- 0.08 / 0.002 


-79.17 


-29.82 / 3.74 


76/ SYM 


5 


0.06 / 0.003 


-80.02 


-28.66 / 3.83 



29 



B. COMPARISON BY PHASE SHIFTER BITSIZE 



In this section, the phase shifter bit size is varied 
while all of the other array parameters are held fixed. For 
all algorithms increasing the number of bits improved the 
sidelobe level, null depth, and beam pointing error. Therefore 
the data presented in this section will be limited to the 
symmetric running sum algorithm because the trends for the 
other roundoff algorithms are similar. The data is summarized 
in Table 8. Figures 13, 14, and 15 show the radiation pattern 

TABLE 8. COMPARISON OF RADIATION CHARACTERISTICS OF ANTENNA 



ARRAYS USING SYMMETRIC RUNNING SUM ALGORITHM 
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# Of ; 
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: ired... 
SLL' 
<dB)" 
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.Variance • ; 
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DeptI); 
Mean' •'kIP 
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i-iLevel 

Variance (dB) 
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for a 76 element array using 3, 4, and 5 bit phase shifters, 
respectfully. 
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Figure 13. 76 Element Array Bayliss Difference Beam for a 3 
Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 14. 76 Element Array Bayliss Difference Beam Radiation 
Pattern for a 4 Bit Phase Shifter using Symmetric Running Sum 
Roundoff 
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Figure 15. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figures 16, 17, and 18 compare the pointing error, null 



depth, and sidelobe level as a function of bitsize. 
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Figure 16. Pointing Error versus Scan Angle for a 76 Element 
Array with 3, 4, and 5 Bit Phase Shifters using Symmetric 
Running Sum Roundoff 
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Figure 17. Null Depth versus Scan Angle for a 76 Element Array 
with 3, 4, and 5 Bit Phase Shifters using Symmetric Running 
Sum Roundoff 
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Figure 18. Maximum Sidelobe Level versus Scan Angle for a 76 
Element Array with 3, 4, and 5 Bit Phase Shifters using 
Symmetric Running Sum Roundoff 
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C. COMPARISON BY THE NUMBER OF ARRAY ELEMENTS 



In this section, the number of array elements is varied 
for each bitsize and roundoff algorithm. Figures 19, 20, and 
21 are plots of 26, 50, and 7 6 element arrays with a 3 bit 
phase shifter using the running sum roundoff algorithm. The 
results are similar to those in the previous section. It can 
be seen in these figures, that as the number of radiating 
elements increased, the sidelobe level increase above the 
desired level is smaller, as expected. A larger sampling of 
the various combinations of the number of radiating elements, 
bitsize of the phase shifter, and type of roundoff algorithm 
is included in Appendix A. Appendix B has plots of null depth, 
pointing error, and sidelobe level as a function of scan 
angle . 

D. FORMULAS FOR EVALUATION OF THE AVERAGE SIDELOBE LEVEL 

It has been shown that the maximum sidelobe level of an 
antenna array is determined by the variance of the 
quantization error [Ref. 3: p. 2]. By taking the expected 
value of the radiation pattern, an expression for the mean 
error level is obtained. This was first done for reflectors by 
Ruze [Ref. 3]. A similar approach for discrete arrays of 
elements yields. 



P .^2 cos (6) 

^norm t - - ^ 



N 



(V-1) 
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Figure 19. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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Figure 20. 50 Element Array Difference Beam Radiation Pattern 

for a 3 Bit Phase Shifter using Running Sum Roundoff 
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Figure 21. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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where is the expected value of the normalized power 

pattern (i.e., radiation pattern) and 

¥ is the variance of the quantization error 
0 is the pattern angle 
N is the number of radiating elements 
is the error free power pattern 
The quantization error variance is calculated using, 



if 2 = 



N 



E 



1^1 



(V-2) 



where is the quantization error at the nth element and if is 
the mean of the quantization error. 

Equation V-1 gives an estimate of the sidelobe level based 
on the variance of the quantization error resulting from a 
roundoff algorithm. Since the statistics of the quantization 
error for all the roundoff algorithms is determined primarily 
by the bitsize, the variances associated with each algorithm 
did not differ by large amounts. Equations V-1 and V-2 
explain why the sidelobe levels were essentially independent 
of roundoff algorithm. Figure 22 is a plot of the average 
sidelobe level predicted by the first term on the right hand 
side of V-1 and a power pattern calculated using weighted 
random roundoff. 
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Figure 22. Plot of Predicted Average Sidelobe Level and Actual 
Power Pattern using Weighted Random Roundoff 
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E. FORMULAS FOR CALCULATING THE BEAM POINTING ERROR 



An accurate determination of target location by a radar is 
made possible by the difference beam null position. A closed 
form expression has been derived by Frank and Ruze [Ref. l:p. 
11-37,38] for the beam pointing error when regular roundoff is 
used. The difference pattern can be written as 



where now includes the phase necessary to scan the beam to 
as well as any quantization error. For a difference 
pattern at 0^, the beam pointing error, 60 is 



If the phases are nearly those required to produce a null at 
0g, the sine can be replaced by its argument, and equation V-4 
reduces to. 




(V-3) 



-t-sin( ”^^ sin (9) -i|> w) ) . 
^ 2 



N 



60=V sin( 

n=l ^ 



sin(0)-i(r^) . 



(V-4) 



N 



60=sin ^ (0) 




(V-5) 
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Figure 23 is a plot of the predicted pointing error using 
equation V-5 versus the calculated pointing errors for regular 
roundoff . 




Figure 23. Predicted Pointing Error (*) and Calculated 
Pointing Error (+) versus Scan Angle for a 50 Element Array 
with a 4 Bit Phase Shifter using Regular Roundoff 



44 



V. CONCLUSIONS 



The running sum roundoff algorithms reduce the beam 
pointing error by a factor of 10 compared to the other 
roundoff algorithms, without degrading the maximum sidelobe 
levels or gain performance. The symmetric running sum roundoff 
algorithm does not seem to offer any significant performance 
advantage over the running sum algorithm. The only possible 
advantage the symmetric algorithm may offer is in beam 
scanning rate. Since the roundoff algorithm is only applied to 
N/2 elements as opposed to N elements, a faster beam scanning 
rate may result. 

Even though the running sum algorithm yielded essentially 
zero beam shift and 60 to 80 decibel null depths, these 
numbers will not be achieved in practice due to manufacturing 
and assembly errors. Phase shifters are designed to a set of 
electrical and mechanical tolerances. In a 4 bit phase 
shifter, for example, the least significant bit may not be 
exactly 22.5®. It could be slightly more or less depending on 
how that individual phase shifter's circuits were fabricated 
and assembled. These errors tend to be random from phase 
shifter to phase shifter and therefore their effects can be 
estimated using statistical approaches such as those described 
in Chapter IV. 
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APPENDIX A - PLOTS OF RADIATION PATTERNS AND QUANTIZATION 
ERRORS 

This section will allow the reader to examine the detailed 
effect of the number of elements, bitsize, and the roundoff 
algorithm on the radiation pattern and quantization error. The 
radiation patterns were arbitrarily chosen for a scan angle of 
5 degrees using 26, 50, and 76 elements with 3, 4, and 5 bit 
phase shifters using regular, weighted random, and running sum 
algorithms. The data is contained in figures 24 through 62. 




Figure 24. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter with Perfect Phase 
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Figure 25. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Regular Roundoff 
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Figure 26. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 27. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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Figure 28. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Symmetrical Running Sum 
Roundoff 
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Figure 29. 26 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Regular Roundoff 
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Figure 30. 26 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 31. 26 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Running Sum Roundoff 
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Figure 32. 26 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 33. 26 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Regular Roundoff 
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Figure 34. 26 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 35. 26 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Running Sum Roundoff 
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Figure 36. 26 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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for a 3 Bit Phase Shifter with Perfect Phase 
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Figure 38. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Regular Roundoff 
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Figure 39. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 40. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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Figure 41. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Symmetrical Running Sum 
Roundoff 
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Figure 42. 50 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Regular Roundoff 
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Figure 43. 50 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 44. 50 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Running Sum Roundoff 
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Figure 45. 50 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 46. 50 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Regular Roundoff 



68 




-100 -80 -60 -40 -20 0 20 40 60 80 100 

DEGREES 



in 

hi 

u 

(T 

0 

llJ 

Q 



q: 

0 

q: 

q: 

hi 

• 

h 

z 

< 

D 

0 




0 5 

PT. ERROR = 0.043 



10 15 20 

ELEMENT NUMBER 



25 30 

MAXSLL= -22.93 
NULL AT -67.47 



Figure 47. 50 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 48. 50 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Running Sum Roundoff 
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Figure 49. 50 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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for a 3 Bit Phase Shifter with Perfect Phase 
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Figure 51. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Regular Roundoff 
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Figure 52. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 53. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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Figure 54. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 55. 76 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Regular Roundoff 
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Figure 56. 76 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 57. 76 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Running Sum Roundoff 
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Figure 58. 76 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Symmetric Running Sum Roundoff 



80 




-100 -80 -60 -40 -20 0 20 40 60 80 100 

DEGREES 




0 0 



10 20 30 40 50 



60 



70 



80 



PT. ERROR = 0 



ELEMENT NUMBER 



MAXSLL= -24.01 
NULL AT -70.53 



Figure 59. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Regular Roundoff 
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Figure 60. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 61. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Running Sum Roundoff 
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Figure 62. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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APPENDIX B 



RADIATION PATTERN CHARACTERISTICS 



This section will allow the reader to examine the 
radiation pattern characteristics for scan angle from 0 to 60 
degrees . The plots of the radiation pattern characteristics 
are grouped by array size, each using the same roundoff 
algorithm. In each figure the performance of 3,4, and 5 bit 
phase shifters is compared, in that order from top to bottom. 
The data is presented in figures 63 through 98. 




Arrays with 3, 4, and 5 Bit Phase Shifters using Regular 

Roundoff 
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Figure 64. Null Depth versus Scan Angle for 26 Element Arrays 
with 3 , 4, and 5 Bit Phase Shifters using Regular Roundoff 
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Figure 65. Maximum Sidelobe Level versus Scan Angle for 26 
Element Arrays with 3 , 4 , and 5 Bit Phase Shifters using 

Regular Roundoff 
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Figure 66. Point Error versus Scan Angle for 26 Element Arrays 
with 3, 4, and 5 Bit Phase Shifter using Weighted Random 

Roundoff 
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Figure 67. Null Depth versus Scan Angle for 26 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Weighted Random 
Roundoff 
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Figure 68. Maximum Sidelobe Level versus Scan Angle for 26 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Weighted Random Roundoff 
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Figure 69. Pointing Error versus Scan Angle for 26 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Running Sum 
Roundoff 
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Figure 70. Null Depth versus Scan Angle for 26 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Running Sum Roundoff 
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Figure 71. Maximum Sidelobe Level versus Scan Angle for 26 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 

Running Sum Roundoff 
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Figure 72. Pointing Error versus Scan Angle for 26 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Symmetric 
Running Sum Roundoff 
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Figure 73. Null Depth versus Scan Angle for 26 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Symmetric Running 
Sum Roundoff 



95 



-10 



CO 

TJ 



(/) 

X 

< 

2 



-20 



-30 




10 



20 30 40 

SCAN ANGLE (DEGREES) 



50 



60 




SCAN ANGLE (DEGREES) 




0 10 20 30 40 50 60 

SCAN ANGLE (DEGREES) 



Figure 74. Maxiinuin Sidelobe Level versus Scan Angle for 26 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 

Symmetric Running Sum Roundoff 
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Figure 75. Pointing Error versus Scan Angle for 50 Element 
Arrays with 3, 4 and 5 Bit Phase Shifters using Regular 

Roundoff 
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Figure 76. Null Depth versus Scan Angle for 50 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Regular Roundoff 
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Figure 77. Maximum Sidelobe Level versus Scan Angle for 50 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Regular Roundoff 
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Figure 78. Pointing Error versus Scan Angle for 50 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Weighted 
Random Roundoff 
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Figure 79. Null Depth versus Scan Angle for 50 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Weighted Random 
Roundoff 
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Figure 80. Maximum Sidelobe Level versus Scan Angle for 50 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 

Weighted Random Roundoff 
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Figure 81. Pointing Error versus Scan Angle for 50 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Running Sum 
Roundoff 
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Figure 82. Null Depth versus Scan Angle for 50 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Running Sum Roundoff 
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Figure 83. Maximum Sidelobe Level versus Scan Angle for 50 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 

Running Sum Roundoff 
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Figure 84. Pointing Error versus Scan Angle for 50 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Symmetric 
Running Sum Roundoff 
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Figure 85. Null Depth versus Scan Angle for 50 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Symmetric Running 
Sun Roundoff 
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Figure 86. Maximum Sidelobe Level versus Scan Angle for 50 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Symmetric Running Sum Roundoff 
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Figure 87. Pointing Error versus Scan Angle for 76 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Regular 

Roundoff 
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Figure 88. Null Depth versus Scan Angle for 76 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Regular Roundoff 
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Figure 89. Maximum Sidelobe Level versus Scan Angle for 76 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Regular Roundoff 
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Figure 90. Pointing Error versus Scan Angle for 76 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Weighted 

Random Roundoff 
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Figure 91. Null Depth versus Scan Angle for 76 Element Arrays 
with 3, 4/ and 5 Bit Phase Shifters using Weighted Random 
Roundoff 



113 



-10 





0 10 20 30 40 50 60 

SCAN ANGLE (DEGREES) 




Figure 92. Maximum Sidelobe Level versus Scan Angle for 76 
Element Arrays with 3, 4 , and 5 Bit Phase Shifters using 

Weighted Random Roundoff 
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Figure 93. Pointing Error versus Scan Angle for 76 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Running Sum 
Roundoff 



115 




0 10 20 30 40 50 60 

SCAN ANGLE (DEGREES) 




0 10 20 30 40 50 60 

SCAN ANGLE (DEGREES) 




Figure 94. Null Depth versus Scan Angle for 76 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Running Sum Roundoff 
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Figure 95. Maximum Sidelobe Level versus Scan Angle for 76 
Element Arrays with 3, A, and 5 Bit Phase Shifters using 
Running Sum Roundoff 
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Figure 96. Pointing Error versus Scan Angle for 76 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Symmetric 
Running Sum Roundoff 
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Figure 97. Null Depth versus Scan Angle for 76 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Symmetric Running 
Sum Roundoff 
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Figure 98. Maximum Sidelobe Level versus Scan Angle for 76 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Symmetric Running Sum Roundoff 
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APPENDIX C 



COMPUTER PROGRAMS 



This section includes all the major programs used in 
compiling the data for this thesis. The programs for plotting 
and saving the data are not included in this appendix. All 
programs were written in MATLAB. 

A. MAIN PROGRAM 

Computes the Radiation pattern of a linear array with either 
a Bayliss, Taylor, or a Uniform Distribution using any size 
of Antenna Array with any Bitsize Phase Shifter. 

begin_scan_angle=0 ; % Initializing the initial 

end_scan_angle=60 ; % conditions for the program 

start_angle=-90 ; 
stop_angle=90 ; 
delta=. 5 ; 

num_elements= [ 26 50]; 
phase_algorithm=l ; 
bit_size=[3 4 5]; 
selection=7 ; 

spacingjper_wavelength=. 3 ; 
theta_step=l ; 

wl=l; w2=l; w3=l; w4=l; w5=0; 
dl=0; d2=l; d3=0; 
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rl=0; r2=l; r3=0; r4=0; 



desired_algorithms= [ 2 ] ; 

sll=25; 

nbar=5 ; 

elem_counter=0 ; 
w6=[ ] ; 
w7= ' no ' ; 

mainmenu % Calls a pulldown menu program for parameters 
% modifications 

clc 

bk=2*pi ; 
r=pi/180 ; 
del=. 5 ; 

it=round ( (stop_angle-start_angle) /delta) +1 ; 
pe=0 ; 

while selection == 6 

al=length (num_elements) ; % The "a" varibles are used to 

a2=length (bit_size) ; % initialize the size of the 

a3=length (desired_algorithms) ; % arrays depending on the 

% selections 

a4=al*a2*a3; % made for the particular 

% simulation 

a5=length (begin_scan_angle : theta_step: end_scan_angle) ; 
nel=max (num_elements) ; 

pterr=zeros (a4 , a5) ; % Point Error Matrix 

pred_err=zeros (a4 , a5) ; % Predicted Pointing Error Matrix 
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nul ldpth=zeros ( a4 , a5 ) ; 
sllmax=zeros ( a4 , a5 ) ; 
ang=zeros ( a5 , it ) ; 
e=zeros ( a5 , it ) ; 
phss=zeros(a5, 16) ; 
amp=zeros ( a5 , nel ) ; 
xsi=zeros(a5,nel) ; 
qph=zeros ( a5 , nel ) ; 
trunc=zeros ( a5 , nel ) ; 
scan=zeros ( 1 , a5 ) ; 
angl=zeros ( a5 , 2001 ) ; 

el=zeros(a5,2001) ; 

count=0; 



% Null Depth Matrix 
% Maximum Sidelobe Level Matrix 
% Angle Matrix 
% E Field Matrix 
% Phase State Matrix 
% Element Amplitude Matrix 
% Phase Matrix 
% Quantized Phase Matrix 
% Quantization Error Matrix 
% Scan Angle Matrix 
% Angle Matrix for -1 to +1 degrees 
% about Scan angle 
% E field Matrix for -1 to +1 
% degrees about Scan Angle 



for element=num_elements 

elem_counter=elem_counter+l ; 
for phasemod=desired_algorithms 
for nbits=bit_size 
count=count+l ; 
ths_count=0; 
plot_counter=l ; 

for ths=begin_scan_angle : theta_step : end_scan_angle 
ths_count=ths_count+l ; 

minl=0; % Initializing Variables to Determine 
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maxl=-120; % Pointing Error, Null Depth, and 

sllmaxl=-120 ; % Maximum Sidelobe Level 

sllmax2=-120; 

scan ( ths_count ) =ths ; 

% BAYLISS Coefficients for Difference Beams 



if dl == 1 

taylor 

elseif d2 == 1 
bayliss 

elseif d3 
uniform2 

end 



% Calling the Appropriate 
% program for the Selected 
% Distribution 



1 



% Determining the maximum Amplitude to Normalize the 
% Amplitudes 

am=max(abs(amp(ths_count, l;num_elements(elem_counter) ) ) ) ; 

if am == 0 
am =1; 
end 

amp( ths_count, : )=amp( ths_count, : )/am; 

% Scan Angles in Degrees 
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us=sin(ths*r) ; 

psis=bk*spacing_per_wavelength*us; 

% Generate exact Phase Required at Each Element. When 
% PHASEMOD =0, the Phase is exact and the Phase Slope is 
% Continuous ( no Modulo 360 or 2pi. Note that Positive 
% Scan Corresponds to Increasing Phase Lag with Increasing 
% n . 

for i=l : num_elements (elem_counter ) 
if sym == 1 & phasemod == 3 

xsi ( ths_count / i )=( -1 ) * ( 2*i-(num_elements (elem_counter )+l ) ) /2 

*psis ; 

else 

xsi ( ths_count ,i) = (-l)*(i-l) *psis ; 

end 

end 

if xsi ( ths_count, 1 ) 

< xsi ( ths_count, num_elements (elem_counter ) ) 
xmin=xsi ( ths_count , 1 ) ; 
else 

xmin=xsi ( ths_count , num_elements ( elem_counter ) ) ; 

end 

% Normalize to the Most Negative Phase so that all xsi(i)> 0 
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xsi(ths_count, : )=xsi(ths_count, : )-xitiin; 
qph( ths_count/ ; )=xsi(ths_count, : ) ; 

% number of BITS is nbit; Number of Phase States is 
% 2'' (nbit); Each Bit is 360/(2^nbit) Degrees 

if phasemod > 0 
ns tates=2 '' ( nbits ) ; 
bitsize=2*pi/nstates; 

for i=l:nstates 

phss (ths_count, i)=bitsize* (i-1 ) ; 
end 

if phasemod == 1 % Calling the Desired 

truncate % Roundoff Algorithm 

elseif phasemod == 2 
roundoff 

elseif phasemod == 3 
runround 

elseif phasemod == 4 
other 

end 

end 

% Determining the Quantization Error 
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trunc(ths_count, : )=qph(ths_count, : )-xsi(ths_count, : ); 

% Begin Pattern Loop 
for i=l:it 

th=start_angle+( i-1 ) *delta; 
ang(ths_count, i)=th; 
u=sin(th*r ) ; 

psi=bk*spacing_per_wavelength*u; 

sum=0; 

% sum to get Array Factor 

for nn=l :num_elements(elem_counter) 

q phase=qph( ths count, nn) ; 

if (nn > (num_elements (elem_counter ) /2 ) ) & (d2 == 1) 

q phase=q phase+pi; 

end 

sum=sum+amp(ths_count,nn) *exp( j*(psi*(nn-l )+q_phase) ) ; 
end 

% Normalize the Pattern to the Number of Elements 
ev=abs ( sum) /num_elements (elem_counter ) ; 
e(ths_count, i)=ev''2; 

e(ths_count , i)=10*logl0(e(ths_count, i) ) ; 
end 

if d2 == 1 % Refining the Radiation 
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end 



nuldepth 

elseif dl ==1 j d3 ==1 
max lobe 



% increment to ,001 Degrees 
% to determine the Point 
% Error, Null Depth, and 
% Maximum Sidelobe Level 



pterr ( count, ths_count)=pe; % Saving the Pointing Error 

nulldpth ( count, ths_count)=minl-maxl; % and Null Depth to the 

% appropriate Matrix 



err_pred 



% Determining the Predicted 
% Pointing Error 



plot_nol 



end 



% Calling the Plotting Program 
% for the Radiation Pattern and 
% Quantization Error 



plot_no0 



mat save 



% Calling the Plotting Program 
% for the predicted versus 
% calculated pointing Error 

% Saving the data to plot the 
% various combinations of 
% pointing error, null depth, and 
% Max. Sidelobe level vs. Scan 
% Angle 
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plot(0: 60,pterr(count, : ) ) , grid, xlabel('SCAN ANGLE (DEG) ') , 
ylabel( 'POINTING ERROR (DEGREES) ') ,title(labell) 
meta scantot 

plot (0:60, nulldpth( count , : ) ) , grid,xlabel ( 'SCAN ANGLE (DEG) ' ) , 
ylabel('NULL DEPTH (dB) ' ) ,title(labell) 
meta scantot 

plot(0: 60,sllmax(count, : ) ) ,grid,xlabel( 'SCAN ANGLE (DEG)'), 

ylabel( 'MAXIMUM SLL (dB )'), title ( labell ) 

meta scantot 

!gpp scantot 

!del scantot. met 

savefile % Saving the data 

end 

end 

end 

if w3 == 1 

plotcomp % Plotting the the data save by “matsave" 

end 

selection =7; 
end 
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B. DETERMINATION OF THE POINTING ERROR 



This program determines the pointing error, null depth, and 
maximum sidelobe level for a Taylor Distribution and 
normalizes the radiation pattern to 0 decibels. 

del=. 5; 

it2=-l* ( start_angle+l ) ; 
it3=-l* ( start_angle-l ) ; 

for i=( ( it2+ths ) /del ) : . 002 : ( ( it3+ths ) /del ) 
th=start_angle+i*del ; 

angl ( ths_count , (i-( ( it2+ths ) /del ) ) *500+1 )=th; 
u=sin(th*r) ; 

psi=bk*spacing_per_wavelength*u; 

sum=0; 



% sum to get Array Factor 
for n=l ; num_elements (elem_counter ) 
q_phase=qph ( ths_count , n ) ; 
if (n > (num_elements (elem_counter ) /2 ) ) 
q_phase= q_p hase+pi; 
end 

sum=sum+amp( ths_count ,n) *exp( j* (psi* (n-1 ) +q phase) ) / 
end 

% Normalize the Pattern to the Number of Elements 
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evl=abs ( sum ) /num_elements ( elem_counter ) ; 
dummy=evl*evl; 

if abs( dummy) < le-9 
dummy=le-9; 
end 



el ( ths_count, ( i-( ( it2+ths ) /del ) ) *500+1 ) =10* loglO (dummy ) ; 

evl=0; 

dummy =0 ; 

% Determining the pointing error and null depth 

if (el(ths_count, (i-( (it2+ths)/del) )*500+l) ) < (mini) 

minl=el ( ths_count , (i-( ( it2+ths ) /del ) )*500+l) ; 

theta=th; 

pe=theta-ths ; 

end 

end 



for i=l:it 

if e(ths_count, i) > maxi 
maxl=e( ths_count, i) ; 
thetamax=ang( ths_count , i) ; 
k=i; 



% Determining the location 
% of the first maximum 
% of the differenece beam 
% radiation pattern 



end 
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end 



m=k; 

max2=maxl ; 

% Determining the location of max null if the maximum of the 
% radiation pattern lies on the positive side of the null 

if thetamax > theta 
for i=m:-l:l 

if e( ths_count, i) <= max2 

max2=e ( ths_count , i ) ; 

kl=i; 

else 

break 

end 

end 

b2=k; 

% Determining the location of the second maximum of the 
% difference beam radiation pattern 

for i=kl:-l:l 

if abs(max2) >= abs(e(ths_count,i) ) 

max2 =e ( ths_count , i ) ; 

k2=i; 



132 



else 



break 

end 

end 

bl=k2; 

end 

% Determining the location of maximum null if the max of the 
% radiation pattern lies on the negative side of the null 

if theta > thetamax 

for i=m:it 

if e( ths_count, i) <= max2 

max2 =e ( t hs_count , i ) ; 

kl=i; 

else 

break 

end 

end 

bl=k; 



% Determining the location of the second maximum of the 
% difference beam radiation pattern 

for i=kl:it 
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if abs(max2) >= abs(e(ths_count, i) ) 



max2=e( ths_ 

k3=i; 

else 

break 

end 

end 

b2=k3; 

end 


count , i ) ; 


amax=0 ; 

for i=(bl-2 ):-!:! 


% Determining the location 


if e(ths_count, i) < amax 


% of the first null on the 


amax=e ( ths_count , i ) ; 


% negative side of the scan 


i=i; 

else 

break 

end 

end 


% angle 


amax=0 ; 

for i=(b2+2 ) : it 


% Determining the location 


if e(ths_count, i) < amax 


% of the first null on the 


amax=e(ths count, i); 


% positive side of the scan 


ll=i; 

else 


% angle 
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end 



break 



end 



for i=l:l 


% Determining the maximum 


if e( ths_count , i) > sllmaxl 


% sidelobe level on the 


sllmaxl=e(ths_count, i) ; 


% negative side of the 


end 

end 


% scan angle 


for i=ll:it 


% Determining the maximum 


if e( ths_count , i) > sllmax2 


% sidelobe level; on the 


sllmax2=e( ths_count, i) ; 


% positive side of the 


end 

end 


% scan angle 



if (abs(amax)+3) >= abs ( sllmax2 ) 
sllmax2=-1000; 

end 



if sllmaxl > sllmax2 


% Determining the maximum 


s 1 Imax ( count , t h s_count ) =s 1 Imax 1 - 
else 


-maxi; % sidelobe level 


s 1 Imax ( count , ths_count ) =s 1 lmax2 ■ 
end 


-maxi; 
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e(ths_count, : )=e(ths_count, ; )-maxl; % Normalizing the 

% Radiation 

yb=min(el (ths_count, : ) ) ; % pattern 

if yb <= -89.90 

[ xa , xb ] =min ( e 1 ( ths_count , : ) ) ; 
el ( count, xb)=-90+maxl ; 
minl=el ( ths_count,xb) ; 

end 

el (ths_count, : )=el (ths_count, : )-raaxl; 

x=( (it2+ths)/del)+l; 
y=( ( it3+ths ) /del )+l ; 

C. BAYLISS DISTRIBUTION 

Computes Linear BAYLISS Coefficients for Difference Beams 
mu=zeros (1,11); 
z=zeros (1,10); 
b=zeros (1,10); 

% Linear Bayliss Distribution 

nbarl=nbar-l ; 
for i=l:ll 
mu(i)=i-.5; 
end 
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sll25=sll-25; 

Z(l)=1.87+sll25* .038; 
z(2)=2.5+sll25*.016; 
z(3)=3.35+sll25* .019; 
z(4)=4 .25+sll25*.016; 
a=l .45+sll25*.042; 

for ns=5:nbar 
z(ns)=sqrt(a*a+ns^2 ) ; 
end 

s igma=mu ( nbar + 1 ) / z ( nbar ) ; 

for iTuns=l:nbar 
bb=l ; 

for ns=l:nbarl 

bb=bb* ( 1- (mu (rams ) / ( sigraa*z ( ns ) ) ) '"2 ) ; 
end 



bbb=l; 

for lls=l:nbar 

if ((lls-mms) < 0) j ((lls-rams) > 0) 
bbb=bbb*( l-(mu(mms)/mu(lls) )^2) ; 
end 



end 
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bes=(-l)^(mins) ; 

b ( mms ) =mu ( mms ) 2 /bes *bb/bbb ; 

end 

for alf=l; (num_elements (elem_counter ) /2 ) 
rho=( 2*alf-l ) /num_elements(elem_counter ) ; 

for lls=l:nbar 
pmu=mu ( 11s ) *pi*rho ; 
bes=sin(pinu) ; 
gg=gg+b( 11s ) *bes ; 
end 

amp ( ths_count , num_e laments ( elem_counter ) /2+1-alf ) =abs ( gg ) ; 
amp ( ths_count , num_elements ( elem_counter ) /2+alf ) =abs ( gg) ; 
end 
end 

D . ROUNDOFF ALGORITHMS 
1. Regular Roundoff 

% Computes the Phase at each Array element by Truncation 

for i=l :num_elements(elem_counter) 

nsteps=round(xsi(ths_count, i)/bitsize) ; 

xs=nsteps*bitsize; 

qph ( ths_count , i ) =xs ; 

end 
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2 . Weighted Regular Roundoff 

% Program Quantizes Phase Slope using Random Roundoff 

for i=l :num_elements(elem_counter) 

nsteps=f ix(xsi ( ths_count , i ) /bit size) ; 

plo=nsteps*bitsize; 

phi=plo+bitsize; 

qph ( ths_count , i ) =plo ; 

xn=rand; 

% Criteria for Weighted Random Roundoff 

pnorm=(xsi (ths_count , i ) -plo) /bit size; 
r nd=pnorm+xn ; 

if (rnd > 1) 

qph ( ths_count , i ) =phi ; 

end 

end 

3. Running Sum Roundoff 

% Computes the Phase at each Array element by "the Running 
Sum Method" and "Symetric Running Sum Method" 

s=0; 

if sym == 1 

nloop=num_elements ( elem_counter ) /2 ; 
else 

nloop=num_elements (elem_counter ) ; 
end 
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for i=l:nloop 

nsteps=fix(xsi(ths_count, i)/bitsize) ; 

phi_lo=nsteps*bitsize; 

phi_hi=phi_lo+bitsize; 

sone=s+phi_lo-xsi ( ths_count , i ) ; 
stwo=s+phi_hi-xsi(ths_count, i) ; 
if abs(stwo) < abs(sone) 
s=stwo; 

qph ( ths_count , i ) =phi_hi ; 
if sym == 1 

qph ( ths_count , num_e laments ( elem_counter ) +l-i ) =-phi_hi ; 

end 

end 

if abs(sone) < abs(stwo) 
s=sone; 

qph ( ths_count , i ) =phi_lo ; 
if sym == 1 

qph ( ths_count , num_elements ( elem_counter ) +l-i ) =-phi_lo ; 

end 

end 

end 
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